Kadanoff-Baym equations and approximate double occupancy in a Hubbard dimer. 
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We use the Kadanoff-Baym equations within the framework of many-body perturbation theory to 
study the double occupancies (fiRfhui) in a Hubbard dimer. The double occupancies are obtained 
from the equations of motion of the single-particle Green's function. Our calculations show that 
the approximate double occupancies can become negative. This is a shortcoming of approximate, 
and yet conserving, many-body self-energies. Among the tested perturbation schemes, the T-matrix 
approximation is the only one providing double occupancies which are always positive. 
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This brief note reports of preliminary results on the 
use of the Kadanoff-Baym Equations (KBE)^ to com- 
pute double occupancies. The case discussed here is that 
of a Hubbard dimer, described by the time-dependent 
Hamiltonian: 

H =y^[-a[ <7 a 2cI - a\ a ai a + w(t)n la } + U ^ ^Rt^Ri) (!) 

a H=l,2 

where U is the on-site and interaction, fiR a — a} R(y an a 
and a =t;4- The term wit) is a time-dependent, spin- 
independent field; in H, all parameters are in units of the 
magnitude of the hopping term (the latter is taken =-1). 
In the, t = 0, initial state the dimer contains two elec- 
trons with opposite spin; this remains true at all times, 
since H has no spin-flip terms. 

The study of this simple system reveals a shortcoming 
of the KBE when they are used together with many-body 
approximations (MBA:s). Following a route based on the 
equations of motion of the single-particle Green's func- 
tion, some conserving MBA:s (namely, the second Born 
approximation, BA, and the GW approximation, GWA) 
may give negative double occupancies (fiR^fim), which 
is clearly unphysical. The results obtained with the T- 
matrix approximation (TMA) are exempt from this prob- 
lem. This is no accident, and in fact we provide an ana- 
lytical proof of why, in the TMA, local double occupan- 
cies are intrinsically positive. 

Kadanoff-Baym equations. - The KBE govern the 
time evolution of the non-equilibrium, two-time, single- 
particle Green's function G(l,2) = -i(T 7 [^(1)^(2)] }, 
where 1 denote single particle space/spin and time labels, 
r\G\t\. Here, T 7 orders the times t\,t<i on the Keldysh 
contour— 7 and the field operators are in the Heisen- 
berg picture; the brackets (} denote averaging over the 
initial state (or thermal equilibrium). Showing explic- 
itly only the time labels (matrix notation/multiplication 
in space and spin indices is adopted), and specializing 
to time t\, we have (idt t — h (ii)) G (tx, t*t) = 6(12) + 
J £ (ti,t) G (t, t%) dt. Here h is the single particle Hamil- 
tonian and £, the kernel of the integral equation, is the 
self energy, which accounts for the interactions and is de- 
scribed within a given MBA, as discussed below. The 
initial state is the correlated ground state (we work at 
temperature, 1//3 — > 0), obtained by solving the Dyson 



equation G = Go + Go£[G]G self consistently, with 
(e — h)Go = 1. In practice we solve these equations using 
an algorithm discussed ini 

The MBA:s we consider are all conserving^, i.e. quan- 
tities such as the total energy, the number of particles 
and momentum are constant during the time evolution. 
The self energy in the BA includes all diagrams up to 
second order. The GWA^ amounts to add up all the 
bubble diagrams which give rise to the screened interac- 
tion, W = U+ UPW, where P(12) = G(12)G(21). In 
this case the self energy is £(12) = G(12)W(12). In the 
TMA£ one constructs the T by adding up all the ladder 
diagrams, T = $ - WT, where $(12) = G(12)G(12). 
The expression of the self energy then becomes £(12) = 
U 2 G(21)T(12). For a detailed account of the different 
MBA:s see e.g£. 

Double occupancy from the KBE. - While in prin- 
ciple G gives access only to expectation values of one- 
body operators, its equations of motion (i.e., the KBE) 
permit to obtain also some quantities which involve ex- 
pectation values of two-body operators, such as, e.g., the 
total energy. This route can be also used for the double 
occupancy dp-±; for on-site (Ur 7^ 0) interactions, 



dn = (fiRtfiRi) = - — 



f £(13)G(31+) d3 



(2) 



RR 



Fig. [T] shows the time-dependent double occupancy on 
site 1 for U = 2 and w (t) = w 9 (t), where Wq = 3. We 
see that the ground state (t = 0) has positive double oc- 
cupancies for all approximations while the BA and the 
GWA become negative after some time. This unphys- 
ical behavior is a drawback that these two MBA:s can 
have. According to Fig[T]thc TMA, among the examined 
MBA:s, is the only one giving always positive correlation 
functions. 

Positiveness of the T-matrix. - All the MBA:s used 
in this work are conserving. This, however, does not 
guarantee that other properties (e.g. spectral features^ 
or response functions) obtained with such MBA:s au- 
tomatically satisfy further basic criteria. In fact, for 
our dimer, density correlations and pair correlations for 
GWA and BA may violate the positiveness condition: 
'(l)^ 1 '(2)^(2)^(1)) > . The pair correlation func- 
tion is a rather sensitive measure of the quality of a MBA. 
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FIG. 1: Double occupancy. 



where the spectral function C RR i is given by 



CW(e) 



A RR ,{e')A RR ,(e-e')de'. 



Thus, C RR >(e) is positive (negative) definite for e < 2/i 
(e > 2fi). The Dyson equation for the T-matrix is 
f(e) = "1(e) - ®(e)UT(e) (here $,T, l> are matrices in 
site indices, and U RR i = U R S RR >). This gives 



and 



(1 + & r U)f^ = <P(1 - f/f a ), 



T(e) - f t( e ) = [1 - f^(e)?7][$(e) - $*(£)][! - J7T a (e)] 



It has been known for a long time that the Random- 
Phase-Approximation (RPA)£ in the ground state gives 
pair correlations which, at metallic densities, are strongly 
negative at short distances^. Later on, negative pair cor- 
relations at short distances were also found within the 
self-consistent GWA4. 

Would the situation improve if the response function 
was calculated self-consistently from changes in £ and 
thereby via the Bethe-Salpeter equation? Such response 
function fulfills all macroscopic conservation laws, but 
not necessarily yields pair correlation functions every- 
where positive. For example, the RPA discussed above 
is the response function in the Kadanoff-Baym sense of 
the conserving Hartrcc approximation but it violates the 
positiveness condition. 

We now turn specifically to the TMA, to show why 
for this approximation the pair correlation is manifestly 
positive. In the ground state, for the TMA, the pair cor- 
relation function is a simple contraction of the T-matrix 
itself, (uRffiRi) = -iT RR (t,t+) = -iT^ R (t,t), (here > 
(<) refers to the electron (hole) part) and its positive- 
ness is a consequence of the positiveness of the T-matrix 
spectral function. The basic building block in the TMA 
is $ RR > = —iG RR /(t — t')G RR '(t — t'). In Fourier space 
we have 



$Rfl'(e) 



C RR .{e')de' 



it] sgn(e' — 2/i) 



(3) 



where we have used the identity (1 — UT)(1 + $U) 
The T-matrix has thus the spectral decomposition 



T(e) 



D(e')de' 



ir\ sgn(e' — 1\x) ' 



1. 



(4) 



where D(e) = [1 - f r (e)U]C(e)[l - UT a (e)}. Conse- 
quently, D(e) is positive definite for e < 2/x, and 



{h R \h Ri ) = -iT RR (t,t + ) 



2/i 



D RR (e')de' >0. (5) 



The above result remains valid i) for a general two-body 
interaction u(r\ — r 2 )S(t 1 — fa), provided we use a sym- 
metrized TMA which includes both direct and exchange 
ladder diagrams; ii) in the presence of an external field 
for t > 0. In this case, the proof is similar to the proof 
that manifestly positive =FiE> give manifestly positive 
=rG? via the KBFjii. 
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